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ABSTRACT 

Context. The energetics in photo-dissociation regions are mainly regulated by the balance between the heating from the photo-electric 
effect acting on dust grains, and the cooling via the copious emission of photons in far-infrared lines. The Orion Bar is a luminous and 
nearby photo-dissociation region, which presents to the observer an ideal edge-on orientation in which to study this energy balance. 
Spatially resolved studies of such a nearby system are essential as they enable us to characterise the physical processes that control 
the energetics of the regions and can serve as templates for distant systems where these processes cannot be disentangled. 
Aims. To characterise the emission of the far-infrared fine-structure lines of [CII] (158yum), [01] (63 and 145/jm), and [Nil] (122yum) 
that trace the gas local conditions, via spatially resolved observations of the Orion Bar. The observed distribution and variation of the 
lines are discussed in relation to the underlying geometry and linked to the energetics associated with the Trapezium stars. 
Methods. Herschel/PACS observations are used to map the spatial distribution of these fine-structure lines across the Bar, with a 
spatial resolution between 4" and 11" and covering a total square area of about 120"xl05". The spatial profile of the emission lines 
are modelled using the radiative transfer code Cloudy. 

Results. The Herschel observations reveal in unprecedented detail the morphology of the Bar. The spatial distribution of the [C II] line 
coincides with that of the [O I] lines. The [Nil] line peaks closer to the ionising star than the other three lines, but with a small region 
of overlap. We can distinguish several knots of enhanced emission within the Bar indicating the presence of an inhomogenous and 
structured medium. The emission profiles cannot be reproduced by a single photo-dissociation region, clearly indicating that, besides 
the Bar, there is a significant contribution from additional photo-dissociation region(s) over the area studied. The combination of both 
the [Nil] and [01] 145yum lines can be used to estimate the [CII] emission and distinguish between its ionised or neutral origin. We 
have calculated how much [CII] emission comes from the neutral and ionised region, and find that at least ~82% originates from the 
photo-dissocciation region. Together, the [CII] 158/jm and [01] 63 and 145yum lines account for ~90% of the power emitted by the 
main cooling lines in the Bar (including CO, H2, etc.), with [O I] 63//m alone accounting for 72% of the total. 

Key words. Infrared: general - HII regions - ISM: individual (Orion Bar) - ISM: lines and bands 

1 . Introduction Due to its proximity (415pc. lMenten et al.Ll2007h and edge- 
on orientation, the Orion Bar has been subject to extensive stud- 
Most of the mass in the Interstellar Medium (ISM) is in neu- ies in the literature, where it is usually adopted as the prototypi- 
}—{ tral form. Photo-dissociation regions (PDRs) are neuti-al re- cal PDR template in the study of high mass star forming regions , 
5^ , gions of the ISM at the interface between the stars and opaque and i s often used to test PDR models (e.g., iPellegrini et all 
cores of molecular clouds, where the heating and chemistry 120091) . It is one of the brightest PDRs wit h an FUV radiation field 
are regulated by the penetrating far ultravi olet (FUV) pho- at the ionisati on fron t of 0= ^1-4x10** Go dTielens & H ollenbacR 
tons of the ionising source(s) (e .g. Hollenbach & T ielensLll999t il985; Marc oni et al.L Il998l) . wher e Gn is the mean interstellar 
iBernard-Salas & TielensL l2005l) . The gas in the PDRs cools radiati on field (1.6x10"^ W m'l iHabingt Il968h . iTleTens et all 
mainly via the emission of far-infrai-ed (FIR) fine structure lines il993) showed how the Bar presented a layered distribution of 
such as [Cii], [Oi], and molecules (e.g. CO, H2). The intensity PolycycHc Aromatic Hydrocarbons (PAHs), H2, and CO. This 
of these lines depends on local conditions in the gas, and when distribution is thought to be the result of an extended gas com- 
compared to PDR models can be used to trace the spatial evolu- ponent seen edge-on, of average gas density between lO'^-lO^ 
tion of the physical conditions of the gas across the PDRs. cm"^. Evidence of clumps have also been sug gested to ex- 

flain high density and temperature tracers (e.g.. 'Burto n et all 
990; vanderWerf, 1996). Recently, Rubin et al. (20Hl) pre- 

Send offprint requests to: J. Bernard-Salas, e-mail: sented mid-IR Spitzer observations that extend beyond the Bar, 

jbernard@ias.u-psud.fr from 2.1 to up to 12'. 1 from the ionising star 6^ On C. They 
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Fig. 1. Overlay of the [OI] 63jum PACS observation on a 
Spitzer/IRAC 8//m image of the Orion Bar A 4x4 raster map 
- 16 overlapping footprints - was performed (see Sect. 2). An 
example of a footprint with its 5x5 spaxels is illustrated on the 
top left of the figure together with their dimensions. The map 
covers the region before and behind the Bar (toward 0' Ori C). 



detect ionised material (Ne^ , IP of 41eV) all the way up to the 
boundary of their observations. They also find evidence for a 
decrease in electron density and an increase of the PDR tracers 
([Si II] 34.8 jum, [Fe II] 26.0 jum, and molecular hydrogen) as 
the dis tance from 0' Ori C increases. In the FIR, Herrmann et aD 
(Il997h mapped the Orion region with a 22" to 55" beam using 
the Kuiper Airborne Observatory (KAO). Their derived column 
densities in the Bar for O" and were in agreement with an 
edge-on geometry. 

In this paper, we present the first maps of the [C II] (158//m), 
[OI] (63 and 145jum), and [Nil] (122;um) lines of the Orion Bar 
region from the recently la unched Herschel Space Observatory 
(HSO, iPilbratt et all l2010h . With its access to FIR and superb 
spatial resolution at these wavelengths, the HSO allows us to 
study the spatial distribution of these lines in unprecedented de- 
tail and so improve our understanding of the ionised and neu- 
tral interface. These observations enable us to trace the spatial 
evolution and excitation conditions of the gas across the illu- 
minated interface of the Bar. This paper is part of a Herschel 
study of the Orion Bar by our group and is complemented 
by Herschel/SPIRE studies of the CO and CH ^ line emission 
dHabart et al.', '2010', '2011t iNavlor et al.L 1201 Ol) . and dust emis- 
sion (Arab et al., 2011). 

The paper is organized as follows. The observations, data 
reduction and line measurements are given in Sect. 2. In the fol- 
lowing section (Sect. 3) the spatial morphology of the lines is 
discussed. In Sect. 4 the distribution and correlation of the lines 
are analysed. Sect. 5 presents a model tailored to reproduce the 
distribution of the lines. Discussions on the origin of the [CII] 
emission and the cooling are presented in Sect. 6 and Sect. 7. 
Finally, the conclusions are summarized in the last section. 

2. Observations and Data Reduction 

The observatio n s wer e taken using the PACS instrument 
(iPoglitsch et al.L l2010l) on board the HSO on 24 February 



2010. These observations are part of the Evolution of the 
Interstellar Medium guaranteed time key project (observa- 
tion ID= 1 342 1 9 1 1 52 ) ^ Abergel et al. , 20 1 0). Four fine-structure 
lines were targeted: [CII] at 158;um, [OI] at 63 and 145//m, and 
[Nil] at 122yum. At these wavelengths, and in increasing order, 
PACS provides a spectral resolution of 3420, 1050, 1170, and 
1270 respectively (v~90-285 km s"'). The observations were 
taken during Science Demonstration Phase (SDP) in the - now 
decommissioned - wavelength switching mode. 

To trace the PDR and H II region interface a 4x4 raster map 
was performed. Each raster position (or footprint) is composed 
by 5x5 spatial pixels referred to as spaxels. For each spaxel the 
line is observed in 16 different spectral scans, each with an up 
and down scan. The configuration at the time of the observa- 
tion is shown in Figure 1 , where the raster map at the epoch of 
observation is overlaid on top of an 8 yum IRAC image of the 
Orion Bar. We note that an additional map in the direction of the 
Trapezium stars and overlapping the current observation is cur- 
rently scheduled for observation and will be presented in a future 
paper Given the brightness of the lines, the minimum exposure 
configuration of one cycle and repetition per line was performed. 
For these maps the best sampling is achieved using Nyquist sam- 
pling, which we adopted for the lines measured in the red chan- 
nel ([CII] at 158yum, [OI] 145jum, and [Nil] at 122jum). This 
consists of a raster point and raster line step of 24". and 22". 
respectively. This results approximately in steps of 2/3 and 1/2 of 
the slit size along both directions. For the [O I] 63//m line located 
in the red channel, Nyquist sampling is achieved with 16". and 
14". 5 raster point and line steps. However, to minimise exposure 
time, and given the strength of the line, the same step size as 
the lines for the red channel was adopted. An position, about 
40' from the center of the map, was taken at the beginning and 
the end of the observation (not shown in the figure). This ob- 
servation is not subtracted from the data as the purpose of the 
wavelength switching mode is to cancel out the background by 
determining a differential line profile. The observation, how- 
ever, confirmed that there was no background contamination to 
any of the lines studied. 

The data were processed using version 6. 1 .0 of the reduction 
and analysis package HIPE (Ott, 2010). From level 0.5 to level 
1 the standard procedures were followed. From this level on the 
cubes were further processed, using proprietary tools, to correct 
for drifting and for flux misalignments between scans. Drifting 
is caused by temperature deviations in the telescope which can 
cause the signal to be modified over time. They can be easily 
identified by comparing the up and down scans. To correct for 
this effect we first passed a median filter in the time domain to 
avoid glitches. The continuum was then fitted using a linear fiQ 
to characterise and remove the drift. The flux misalignment be- 
tween the spectral scans can be due to improper dark subtraction 
or flatfield correction. The former will have an additive effect, 
and the latter a multiplicative one. In our case, the scans could 
be well aligned using an additive value which indicates that in 
our observations the cause of the flux misalignments is the dark 
subtraction. We note that the lines are very strong and these cor- 
rections were minor 

At this point the cubes are exported and we used the ID L- 
based software PACSman dLebouteiller et al.L 1201 1. in preph to 
measure the line fluxes (by fitting a Gaussian) and create the fi- 
nal map. In this code, the line fluxes are measured for all the 
spaxels independently. The lines are not spectrally resolved. To 



' In our data the signal was either linear or slightly curved. A poly- 
nomial fit of second degree produced very similar results. 
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create the final map, PACSman recreates an oversampled pixe- 
lated grid of the observations (with 3" pixel resolution) and cal- 
culates the average fractional contribution of the given spaxels 
to the relevant position. The code calculates the statistical uncer- 
tainties on the fly, which include the dispersion in the reduction 
process and the RMS of the fit. These uncertainties are small and 
usually amount to less than 4%, 1%, 3% for the [CII] and [OI] 
lines. For the weaker [Nil] 122yum line, these are higher and 
oscillate between 5-13%. The relative accuracy between spax- 
els given in the manual is 1096, and for the remaining of the 
paper we have adopted the higher of this value or the uncer- 
tainty given by PACSman. We note that absolute flux calibration 
is quoted to be 30% of the peak-to-peak accuracy (scatter around 
the expected flux densities). The final maps for each of the four 
lines are shown in Figure 2. We have c ompared our [C 1 1] and 
[OI] line fluxes with those measured bv lHerrmann et al] (Il997h 
with the Fabry-Perot interferometer FIFI and integrated over a 
55" beam. Their measurements are in good agreement (within 
~30%) with those reported in this paper when convolved to their 
larger beam size. 

3. Spatial Distribution 

The observed spatial distribution of the lines is shown in 
Figure 2, where the ionising star 6*' Ori C is illuminating the 
region from the top right, outside the maps (see also Fig. 4a for a 
sketch of the different regions). The beams sizes are represented 
in the top right of each map and correspond, in increasing order 
of wavelength, to 4". 5, 8". 8, 10", and 11". At the distance of the 
Bar (415 pc), 10" corresponds to a physical scale of 0.02 pc. The 
Bafl is resolved in all four lines and we detect emission all over 
the region probed, including detection of [N n] (ionised gas) be- 
hind the Bar. The range of intensities is not large, with variations 
by factors 4 ([C II] and [N II]) to 8 ([OI]) from the peak (Bar) to 
the fainter (outer) regions. 

The overall morphology of the region can be best discussed 
from the [O I] 63/im map because its smaller point spread func- 
tion (PSF) provides the highest detail (Figure 2). Even if this 
line is more affected by self-absorption, the gain in resolution 
bests this caveat. The emission peaks at about 123" SE of 0' 
Ori C. The most striking feature is the presence of several knots 
of enhanced emission resulting in small scale structures, which 
suggests a clumpy distribution within the Bar These knots, three 
in the south-west and the two in the north-east, are bridged by 
weaker emission, ~16% lower in flux than the knots. In the 
[OI] 63/im image, the smaller knots have diameters between 
~6" and ~10" which, at the distance of the Orion nebula, cor- 
respond to 0.01-0.02 pc. The smaller ones (6") are marginally 
resolved and could result from the superposition of even smaller 
clumps. These knots could be clumps bei ng photo-evaporated 
by the intense FUV penetrating the PDR. iGorti & HollenbachI 
(l2002h studied the effects of the radiation field on the evolu- 
tion of clumps in PDRs. Using their relation between the size 
of the clumps and the column density of the FUV heated region 
(No=2x 10^' cm"^, their Eq. 34), we derive a density of the gas at 

^ From the PACS spectroscopy performance 
and calibration manual. This can be found at 
|http://herschel.esac.esa.int/twiki/bin/view/Public/PacsCalibrationWeb 

^ For the remainder of the paper, the Bar is defined as the region 
where the emission is higher than 75% of the peak emission for carbon. 
This threshold falls between the third and fourth highest contours in 
Figure 3 and is represented by the solid line in Figure 4a. Similarly, 
we use the same definition for the discussion in Sect. 6 about the [N II] 
peak. 



the base of the photo-evaporating flow of 6-10x10'' cm"^. This 
value is consistent with the inter-clump density we use in Sect. 5 
for the modelling. Knowing that the density inside the clumps is 
higher than this, and assuming that the observed size is that of an 
isolated clump, we estimate that the photo-evaporation timescale 
of these clumps should be higher than ~4 500-7 500 yr, depend- 
ing on size. 

This morphology is somewhat mimicked in the [O I] 145//m 
and [CII] 158;um maps, where the lower resolution has washed 
out the clumps of emission. As we move away from the Bar 
(both in front and behind) the emission decreases gradually. 
There is, however, a confined region of brighter emission in the 
north of the map (middle top of the [OI] 63/im map), proba- 
bly due to an increase in density in the PDR which lies behind 
the HII region. Alternatively this could be the result of a geo- 
metrical effect where the inclination of the background PDR is 
steeper. This region is also revealed in the other three maps, al- 
beit with less contrast. The [OI] 63|/m map also points to an 
increased emission at the western edge of the Bar and extending 
north. This excess was also detected in CO emission from the 
ground (e.g.. lLis et al], ll998l). 

The morphology of the [C II] line follows very well that of 
the [O I] lines (especially that of [O I] 145//m) with the Bar peak- 
ing at the same position in the three maps. This is better seen in 
Figures 3c and 3d where the images have been convolved to the 
same resolution as the [CII]158yL(m map (largest PSF), and are 
shown with over-plotted contours. In Figure 3a we also see that 
the PAH emission, as traced by the 8;um IRAC band, is close to 
the [C II] peak but slightly shifted towards the ionising star 0' Ori 
C. The peak of [C ii] emission in the top-left corner of the Bar is, 
however, not followed as well by the PAH emission. In Figure 3d 
we can see that the different intensity levels in the 145/im map 
delineate the Bar very well, whereas in the 63/im map the spread 
in the levels is broader This is probably the result of the higher 
optical depth on the 63/im line, where changes in column density 
have a less pronounced effect than for the 145/im line. 

The peak of the [N II] emission is displaced with respect to 
the other lines by about 12" towards to 0' Ori C (Fig. 2 and 
Fig. 3b). This is expected because with an ionisation potential 
of 14.5 eV the [Nil] line comes from the ionised phase and is 
representative of the ionisation front. There is a small region of 
overlap between the [N II] and [C II] lines but this could be due 
to an orientation effect (Bar slightly tilted). We even note that 
the north-east tip (top-left) of the Bar that is seen in the [C II] 
158/im and [O I] 145//m maps (and also 8//m IRAC map. Figs. 1 
and 2), delineates the boundary of emission of the [N II] line in 
that region. On the western side of the Bar, the [N II] emission is 
somewhat delineated by the above-mentioned emission extend- 
ing north in the [O I] 63/(m map and also detected in CO and 
dust emission. Finally, the gap in the center of the Bar, seen in 
the [O I] and [C II] lines, coincides with the maximum emission 
of the [Nil] llljim line in front of it (Fig. 3b). This could indi- 
cate that the radiation field in this specific spot is higher (or more 
intense) and has been able to photo-ionise the gas. Alternatively 
this could be due to variations in the column density. 

4. Line Correlations and Ratios 

For the purpose of the discussion we have divided the map into 
several regions which are represented in Figure 4a and annotated 
in panel c of the same figure. Summarizing: data points in front 
of and behind the Bar are colour-coded in blue and red, respec- 
tively; plus symbols indicate points that trace the Bar in [C II] 
(those inside the solid contour line), and dots points outside the 
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Fig.2. PACS contour maps of the [CII] (158/zm), [OI] (63 and 145yum), and [Nil] (122//m) lines in units of 10"" W nT^ sr 



In increasing order of wavelength the beams sizes are 4". 5, 8' 
emphasize the detail. 



.8, 10", and 11". Note the offset from zero in the intensity scale to 



Bar in [CII]. In addition, selected points along an adopted cut 
(see §5) have been labeled (orange) and include error bars. 

Figure 4b plots the [Nil] versus the [CII] line fluxes. The 
points before and behind the Bar (blue and red respe ctively) oc- 
cupy diff'erent regions but there are no clear trends. lAbel et all 
( 120051) hinted, from their models, at a correlation between these 
two lines in the ionised medium. The intricate behaviour in our 
figure does not suggest that. On the other hand. Figure 4c il- 
lustrates that the [CII] 158/vm line follows very well the [OI] 
145;um line. A linear fit gives a slope of 1.059+0.005. Variation 
of the [OI] 145;um line is dominated by the column density of 
the gas, and so this correlation indicates that most of the [C II] 
comes from the PDR and not the HII region (see Sect. 6). 

The [OI] 63jum and [CII] 158yL(m lines show also a good 
correlation (Fig. 4d) but it seems that the correlation outside 
the Bar(CII) differs with respect to the other regions. This is 
also seen in Figure 4e where the two oxygen lines are plotted 
against each other. The points behind the Bar (red dots) are dis- 
placed compared to those before the Bar (blue dots), and there 
is a change of slope between points in the Bar (plus signs) and 



outside the Bar (dots). The excitation conditions of these lines 
are similar and therefore this behaviour could reflect the effect 
of opacities that more strongly affect the [OI]63/jm line (as it 
becomes self-absorbed at rather low column densities). 

In the following panels we plot the ratio of several lines 
which, in an optically thin environment, allows us to remove the 
effect of column density in the trends: 

• Figure 4f shows that the [CII](158//m)/[NII](122jum) ratio 
is different before and behind the Bar A map of this line ra- 
tio is illustrated Figure 5b where we clearly distinguish two 
regions, with the points before the Bar, which are more rep- 
resentative of the ionised gas, having a lower ratio. We find 
that a threshold of 15 for this ratio could be used to distin- 
guish between emission from the neutral and ionised region. 

• The [C II] 158//m/[OI] lASfim ratio in Figure 4g decreases as 
we approach the Bar. The trend in the figure of the points 
outside the Bar (dots) can be approximated by a polyno- 
mial fit. This behaviour is again evident in the map ratio 
shown in Figure 5a where points before and behind the Bar 
show the similar values. This suggests similar physical con- 
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been convolved to the PSF at 158yum and are in flux units of 10"^ W nT^ sr"' (except for the IRAC image which is in 10-''MJy sr"'. 
Here the intensity scales start at zero. 



ditions in these two regions because the excitation param- 
eters for both lines are different. This is not expected for 
a single PDR and suggests contamination of a background 
PDR (see Sect. 5) Another important point is that within 
the Bar the ratio is very homogeneous (within 16%) but it 
differs from the regions outside the Bar. We can use the 
[C II] 158jum/[0 1] 145jum ratio to put a constraint on the den- 
sities b y comparing with the predicted va l ues fr om PDR 
models ( Ferland et allll998l:lKaufman et alill999h . We find 
densities that vary from 10^-10^cm~^ in the Bar, and 10-^- 
lO'^cm"-' outside it. 

In the last panel Figure 4h, it is clear that the line ratio 
[O I]63;um/[0 1] 145;um splits before and behind the Bar (be- 
ing lower before the Bar). Given that the excitation condi- 
tions in these regions must be similar (Fig. 4g), we conclude 
that the observed difference is the result of [O I] 63/im be- 
coming self-absorbed; in the Bar, as it is expected, the ratio 
is smaller because the column density is higher than out- 
side the Bar. We note that the gas in the region where the 
[O I]63;um line is optically thick must still be warm to excite 
the [OI]145yL(m line (E„ =326 K). From Figure 5c we see 
that the ratio is indeed not homogeneous within the Bar (as 
it was for the [C II] and [O I] 145/zm, Fig. 5a), indicating that 
opacity effects are present. 



5. Modelling 

We have modelled the line emiss ion across the Bar using 
the radiative transfer code Cloudy (iFerland et al.L Il998l) . This 
code computes the chemistry and radiative transfer at the sur- 
face of a molecular cloud (in our case assumed to be a 
plane-parallel semi-infinite slab) illuminated by FUV photons 
from the ionising star(s). We adopt a Kurucz model for the 
star at 39 600 K (Pellegrini et al .L l2009l) . and ISM abundances 
(Savage &Sembach, 1996; Mev er etal.L [l997, 1998). In order 
to reproduce the spatial stratification we adjust the starting point 
of the HII region and the density. Namely, the ionised region is 
adjusted to start at 0.134 pc (67") from the star. The adopted gas 
density is set to a constant value of 3200 cm dPellegrini et al.L 
2009) in the ionised region, and is t hen coupled with a profile 
density for the PDR as described in lArab et al] (1201 lb . with a 
density scaled to 6x10"^ cm at the peak of the [OI] and [CII] 
and emission (projected distance of 0.246 pc). No tilting is as- 
sumed for the Bar and the depth of the PDR along the line of 
sight is adjusted to 0.35 pc. 

To illustrate the comparison between the observed fluxes and 
the model (convolved to the 158/im PSF) we assume the cut il- 
lustrated in Figure 4a. This cut is also adopted for the discus- 
sion m ou r accom pany mg papers by iHabart et all (l201ll) and 
lArabetal.1 (120111) . It is made to pass through the exciting star 
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Fig. 4. Intensity plots for different combination of lines fluxes (in 10"^ W m"^ sr ') and ratios from the convolved maps. The 
different regions are indicated in the first panel and are labelled in the third, where blue represents points in front of the Bar and red 
behind the Bar Plus symbols indicate points within the Bar(CII). The square data point is the position where the cooling is calculated 
(see Sect. 7). The reference points for the adopted cut (Sect. 5) are plotted in circles and include error bars to give an indication of 
the uncertainties at different regions. In panels c and g the solid lines represent respectively a linear fit (y = -2.28 + 1.06x), and a 
polynomial fit for the fainter region outside the Bar (y = 6.37 - 1 .53x + 0. 12jt:^). 
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Fig. 5. Map ratios for different combination of lines as indicated in the title of each panel. 



6^ Ori C (5h 35m 16.46s, -5d 23m 23. 17s) and an arbitrary point 
in the Bar (5h 35m 21.82s, -5d 24m 59.18s). The chosen point in 
the Bar minimizes effects of other conta minating stars in th e di- 
rection of the cut that can affect the dust (Arab et 3111201 ih . and 
allows us to sample different conditions in the region. In Figure 6 
the solid line represents the output of the model, and the dotted 
line the model plus an offset (see point 3.) From the comparison 
in Figure 6 we find that: 

1 . This simple model does a very good job in reproducing the 
observed spatial stratification of the peak positions of the 
[Nil] line at 111.5", and the [CII] and [OI] lines at 123.5" 
(top and middle panels). This stratification is the result of the 



attenuation of the incident radiation field across the Bar and 
is sensitive to the density profile. 

2. The model also shows a gradual increase of the [Nil] line as 
we approach the Bar (better seen in the bottom panel) which 
is also seen in the observations, while in contrast the oxy- 
gen lines show a more abrupt increase. This is the result of 
the ionisation structure. The model also shows a gradual in- 
crease for the [C II] line but this is not clear from the observa- 
tions. This could be the result of it being more easily excited 
(low temperature and density) compared to the oxygen lines 
and/or contribution of the [C II] line from the ionised phase. 

3. It is striking that the observed profiles do not fall to zero 
as the model (solid lines) does both before and after the 
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Fig. 6. Observed (plus symbols) and modelled (solid and dotted 
lines) profiles along the cut shown in Figure 4 (see Sect. 5) mea- 
sured with respect to the distance to the exciting star 0' Ori C. 
Observed and modelled fluxes are convolved to the 158yum PSF, 
where the dotted line includes an offset due to background emis- 
sion (see text). The first four panels compare for each line the ob- 
served convolved profiles with a Cloudy model, where this has 
been scaled to the peak of the [O I] (145 jjm) emission. The bot- 
tom panel compares the normalized model profiles. In this panel, 
the solid dots represent the position of the reference points in the 
adopted cut (Fig. 4a). 



Bar. This background emission is significant, amounting to 
about ~18% of the peak emission for the [OI] 145//m line 
and ~38% for the other three lines. These offsets have been 
added to the model emission as the dotted lines in the Figure. 
This suggests the presence of another PDR(s), probably from 
the background cavity in the molecular cloud, which must 
produce the base of emission that we measure. This PDR 
is probably not entirely face-on because the strong intensity 
observed, which may result from a limb brightening effect, 
cannot be reproduced by a face-on PDR model. 

4. The model line profiles are narrower than the observa- 
tions. This could be due to the Bar being tilted to the ob- 
server (^j^2}j_an_effec£_already inferred in other studies 
(e.g., iPellegrini et al.L l2009h . Alternatively, a highly struc- 
tured medium would allow FUV radiation to permeate the 
region and heat gas on larger spatial scales. 

5. In order to reproduce the absolute flux, and assuming that the 
base emission from the additional PDR(s) has an additive 
effect on top of the emission of the Bar (dotted line), we 
need to assume a depth of the PDR in the model of 0.3 pc. 
The depth given above is larger than, but comparable to, the 
width of the Bar (~0.06 pc) as expected for a tilted edge-on 
Bar. 




Fig. 7. Line intensity ratios for different combinations of the 
[Nil] (122 and 205yum) and [CII] 158yum lines. Theoretical 
curves are shown by the solid, dotted, and dashed lines. Grey ar- 
eas indicate the range of observed ratios, and the enclosing boxes 
represent the range of uncertainties in the ratios. The square 
point is the ratio of the [Nil] 122/205 fines. 



6. In this model we take into account optical depth effects 
across the PDR0 but not along the line of sight. However, 
we can see that the model (with offset) reproduces the peak 
of emission of the [Nil] 122/im, [OI] 145/im fairly well, 
and the [CII] 158yL(m within factor 1.2. This indicates that, 
except for the [O I] 63//m line, the emission of these lines 
are not dominated by optical effects. 

A more detailed modelling of the region including additional 
PDR/s (even the effects of dumpiness and tilting) and dust emis- 
sion is outside our scope here but will be presented in a future 
paper. Still, it is encouraging to see how a simple model (single 
edge-on PDR) can be used to at least reproduce with a fair degree 
of success the observations of such a complex environment. 



6. [CII] from the PDR and HII region 

In the ISM, the [CII]158;um line is important in the study of 
the cooling and chemistry of PDRs. In extragalactic studies it 
is important for redshift determinations, and for the extent to 
which its luminosity is a measure of the star formation rate (SFR , 



iityj 

IStacev et al.Lll99UlMeiierink et alil2007HLuhman et alil2003h 
With a low ionisation potential of 1 1 .3 eV, this line can originate 
both in the PDR and in the HII region. It is thus important to 
characterise in detail its contribution in different environments. 
Because of the high radiation field and high density in the Orion 
nebula, one expects that most of the [C II] should originate from 
the PDR. The PACS maps show in great detail the distribution 
of the [C II] and [N II] lines where we can investigate this issue. 

We have therefore quantified the contribution of [CII] 
158^(m and [Nil] 122yum to the Bar as traced by the [C II] emis- 
sion (Bar(CII)), and as traced by the [Nil] emission (Bar(NII)). 
These are represented as areas enclosed by the solid and dashed 
lines, respectively, in Figure 4a. This combined region, which 



We assume both thermal and micro-turbulence line width with a 
turbulence velocity of 3 km s"' . 
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we name Bar(CII+NII), encompasses the peak emission of the 
PDR and ionised region. We find that 76% of the [CII] emis- 
sion is coming from the PDR region (Bar(CII)) relative to the 
total emission in the Bar(CII+NII). The remaining 31% comes 
from the Bar(NII) regioifl This 3 1 % conti'ibution to the Bar(NII) 
region could include emission from the ionised region or from 
the additional PDR (see Sect. 5). To estimate this, and ignor- 
ing confusion from tilting effects, we can use the close rela- 
tion of the [CII] line with the [OI]145;um line. For a broad 
range of densities between 10^ and 10^ (typical of PDRs), the 
[C ni/rO II 145/im rati o varies between 0.5 and 2 jFerland et al.L 
ll998VKaufm an et al.L [19991: iLe Petit et"all 1200 6 '). The observed 
[OI]145yL/m in the Bar(NII) is about 2x10"^ W m ^ sr from 
which we predict a PDR [CII] emission between 1 to 4x10"^ 
W m"^ sr"'. We can then compare this range to our observed 
[C II] emission of 4x10"^ W m"^ sr ' in this region. Thus, from 
the original 3 1 % of [C II] in this region, we deduce that either 
all could come from the additional PDR, or the HII region could 
contribute to at most 24%. 

Following the analysis by lOberst et al.l (1201 ll) in the Carina 
nebula, we have compared the expected [C II]/[N II] ratio in the 
ionised medium with the observed values in the Bar(NII) re- 
gion. This can be used to get an additional estimate of the con- 
tribution from the ionised region to the [CII] line. To do this 
we have complemented our observa tions with t he SPIRE-FTS 
data of the [N II]205;um line (HabartetaTL 1201 ih . and have cal- 
culated the theoretical curves for the [Nn]122yum/[NII]205//m, 
[CII]158A<m/[NII]205|/m, an d [C II] 158A<m/rN II] 122//m ratios 
at a temperature of 9000 K (iNieva & Si^mon-Diazl 1201 ih . The 
[CII]/[NII] line ratios depend on the relative ionic abun- 
dance of CII and NIL We use the carbon and nitrogen abun- 
dances adopted in Sect. 5, coupled with the ionisation frac- 
tion of CII and Nil fo r the HII conditions in the Orion neb- 
ula as given by RubinI ([1981 . to derive a relative ionic abun- 
dance of [C II/NII]~1.6. These theoretical curves are shown 
in Figure ^ The [NII]122//m/[NII]205/im ratio can be used 
to determine the electron density. The observed ratio (square 
in the figure) falls in the non-linear regime of the theoreti- 
cal curve (solid line) but can serve to obtain an upper limit 
in the density of 100 cm~^, consis t ent w ith our adopted 
value of 3200 cm"^ (iPellegrini et al.L l2009l) . The observed 
[CII]158/zm/[NII]205//m, and [CII]158/zm/[NII]122//m line 
ratios, grey areas in the figure, are clearly above the theoreti- 
cal curves (dotted and dashed lines respectively). Taking into 
account the uncertainties in the line fluxes (enclosing boxes in 
the figure), the observed ratios indicate that no more than 9% 
of the [C II] 158;um comes from the ionised region. Allowing for 
an uncertainty of about a factor 2 in the relative abundance of 
C II/N II and collisional strengths used to calculate the theoreti- 
cal curves, we can safely conclude that in the Bar less than 18% 
of the [C II] emission originates in the ionised medium. This is a 
slightly lower threshold than our previous estimate. 

lAbeTeFaLl (j2005j) made some modifications to the Cloudy 
code to, among others things, estimate the contribution of sev- 
eral PDR lines from the H II region. They apply their model to 
the starburst galaxy NGC253 where they find that about 30% 
of the [CII] comes fro m the ionised medium. More recently, 
iMookeriea et al.l (1201 ll) using PACS data have mapped the FIR 



' These quantities do not add to 100% because there is an overlapping 
region as it is explained in Sect. 3. 

^ The [CII]158yum/[NII] (122, 205yum) curves can be applied to any 
other object by simply scaling them to the ratio between the relative 
abundance of [C II/N II] in the object and that assumed here (1.6). 



emission of the [C II] 158;um line in spiral galaxy M33. They find 
that between 20 to 30% of this emission comes from the ionised 
medium. Their values are slightly higher than our determina- 
tions, which is not surprising since we are probing dense PDR in 
the Orion Bar while in these galaxies lower density HII regions 
could dominate the emission. We note, however, that this com- 
parison is hampered by the diff'erence between the larger scales 
(where diff'erent regions are mixed) that these studies probe, and 
the small scales probed in our study. 

7. Cooling 

In order to establish the role of the [CII] and [OI] lines in 
the cooling of the PDR we have compared their contribution 
to that of other relevant cooling lines in the mid- and FIR. We 
have included in this calculation the H2 rotational lines from 
the ISO/SWS, and ^^CO rotational lines (J=4-3 to J=21-20), 
■^CO (J=5-4 to J^1 4-13), H2O, and CH+ from the PACS and 
SPIRE instruments (iHabart et al.Ll2010l:lJobhn et alil201 ll) . To 
make this comparison we use the same position (5h35m21s, - 
5° 25' 18") of the H2 measurements by ISO to derive the con- 
tribution of each line (or cascade of lines). This position is in- 
dicated by the white square in Figure 4a. For this comparison 
we convolve all the data to the largest beam size (SPIRE 40")- 
We find that in this region the [C II] and [O I] lines studied in 
this paper contribute 90% of the total power emitted by all these 
lines, with the [OI] 63;um line contributing 72% of this emis- 
sion. The CO, H2, and CH^ contribute, respectively, 5, 4, and 
less than 1%. These estimates accentuates the importance of the 
[O I] 63/im line to the cooling budget from the gas lines in these 
regions. This fact is also highlighted by the relative strength of 
this line compared to the other lines in the region (see Fig. 6). 

8. Summary and Conclusions 

We have presented the first HSO observations of the 
[CII]158/im, [OI]63//m and 145A<m, and [NII]122/im lines of 
the Orion Bar. Its angular resolution has allowed us to map the 
spatial distribution of these lines in unprecedented detail. 

The [CII] and [OI] maps peak at the same position and 
fall close to the peak emission of PAHs (as traced by the 
Spitzer/IRAC 8yL/m band). The [Nil] peaks sHghtly closer to 0' 
Ori C with a small region of overlap with respect to the other 
PDR lines. Within the Bar we can distinguish between knots 
(clumps) of emission, about 0.01-0.02 pc in size (6 to 10"), and 
which are 16% higher in flux than the interclump medium. These 
clumps could be photo-evaporated by the FUV inside the PDR. 
These knots of emission are better seen in the [O I]63/im map as 
this line offers the best resolution (having the smallest PSF), it is 
also seen in the [O I] 145yum and [C II] 158//m maps. 

The [CII]158yL(m correlates very well with the [OI] 145//m 
emission. The [C II] 158;um line does not correlate with the [N II] 
in the ionised region as some studies have suggested. The com- 
bined information on the [N II] and [O I] lines provides a great 
diagnostic to estimate the emission from the [CII]158/im line, 
and to distinguish between an origin in an ionised or neutral re- 
gion. The ratio between the [OI] 145/^m/63/^m lines show the 
effect of the opacities, where the [O I] 63;um line becomes self- 
absorbed at high column densities. 

We have modelled the emission of the lines with the photo- 
ionisation code Cloudy and reproduce the relative position of 
the lines. The emission profiles reveal that, in addition to the 
Bar, there is a significant background emission all over the re- 
gion (present for all four lines). This points to the presence of 
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additional PDR(s). This should be follow up with detail radia- 
tive transfer models to infer the physical conditions and coupled 
to what has been learned from the emission of dust. 

The [CII] line can come from the neutral and ionised 
medium. We have made different estimations of its contribution 
and find that most of the [CII] emission originates in the PDR(s) 
(>82%). Using ancillary ISO and Herschel data we have calcu- 
lated the total power emitted by the most relevant cooling lines 
from the mid- to FIR (fOll, [CII], CO, H2O, CH+). We show 
that the power emitted by the atomic [Cll] 158yum, and [OI] 
63fim and 145jum lines account for 90% of the power emitted in 
the region by all of the cooling lines considered, with the [OI] 
63yum line contributing 72% of the total. This emphasizes the 
predominant role of the latter in the cooling process from emis- 
sion lines of these regions. 
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